Survival of a Diffusing Particle in a Transverse Flow Field 
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We consider a particle diffusing in the y-direction, dy/dt — n{i) where n(t) is Gaussian white 
noise, and subject to a transverse flow field in the ^-direction, dx/dt = f(y), where x > and x = 
is an absorbing boundary. We discuss the time-dependence of the survival probability of the particle 
for a class of functions f(y) that are positive in some regions of space and negative in others. 



PACS numbers: 02.50.-r, 05.40.-a 
I. INTRODUCTION 

In a recent paper [1] we considered a class of stochastic 
processes denned by the equations 



y = v(t) > 
x = f(y), 



(i) 

(2) 



where rj(t) is Gaussian white noise with mean zero and 
correlator (r)(t)r)(t')) = 2DS(t-t'). Equations (1) and (2) 
represent a particle which diffuses (performs a random 
walk) in the y-direction but is subject to a deterministic 
drift in the x-direction with a y-dependent drift veloc- 
ity. We consider the case where there is an absorbing 
boundary at x = 0, and we are interested in the sur- 
vival probability of the particle. Two previously studied 
models are the cases /(y) = y and f(y) = sgn(y). The 
former is equivalent to the random acceleration process, 
x = rj(t), while the latter is the 'windy cliff model in- 
troduced by Redner and Krapivsky [2]. For both these 
models, the survival probability of the particle decays as 
t -1 / 4 for large time t, and it has been argued [1,2] that 
this is a generic result for odd functions f(y). It should 
be noted that, in general, such problems are highly non- 
trivial and exact solutions are available in only a small 
number of cases. 

In ref. [1], we raised the question of how the sur- 
vival probability decays when the function /(y) is not 
an odd function. We considered a class of models where 
f(y) = V +V a for 2/ > 0, and f(y) = -v_(-y) a for y < 0. 
Here the drift takes the particle away from the absorbing 
boundary for y > 0, and towards it for y < 0, but not in 
a symmetrical way: the function f(y) is odd only when 
v + = V- . Using an extension of a technique developed by 
Burkhardt [3] in connection with the random acceleration 
process, we showed that the survival probability decays 
as t~ e for large t, where the exponent 9 (the 'persistence 
exponent') is nontrivial and given by 



27T/3 



■ tan 



■ tan 



~2~ 



(3) 



where (3 = l/(2+cv). This reduces to 9 = 1/4 for w + = u_ 
but 9 lies in the range 0<#<l/2in general. When 



v + > the drift away from the boundary is weakly 
dominant. The decay of the survival probability still has 
a power law form, but 9 < 1/4, while the converse occurs 
for v+ < V-: the drift toward the boundary is weakly 
dominant and 9 > 1/4. The fact that 1/2 is an upper 
bound for 9 is clear when one recalls that the probability 
that the particle stays in the region y > (and therefore 
never encounters the drift towards the absorbing bound- 
ary) decays as t~ x / 2 [4,5]. 

It is worth emphasising that the process defined by 
Eqs. (1) and (2) is non-Gaussian except when f(y) is lin- 
ear. Eq. (3) is a rare example of an exactly calculable 
persistence exponent for a non-Gaussian process. 

In [1] we speculated about the behavior of the system 
when the function /(y) is described by different expo- 
nents, a+ and a_, for y > and y < respectively. 
We argued that for a+ > a_ the drift away from the 
boundary is strongly dominant, leading to a non-zero sur- 
vival probability at infinite time (i.e. 9 = 0), while for 
a + < a_ the drift towards the boundary is strongly dom- 
inant and 9 obtains its maximum value of 1/2. 

In the present work we obtain some exact result per- 
tinent to the latter question. Looking at the case a + > 
ev_, we show explicitly that the survival probability ap- 
proaches a non-zero value at infinite time and we obtain 
a closed-form expression for this value. 

In the second part of the paper we revisit the issue 
of whether the decay exponent 9 takes the value 1/4 for 
all odd functions f{y), as has been assumed up to now. 
We consider models where the flow field takes a periodic, 
banded form with f(y) taking the values v + and — 1;_ in 
alternate bands of width b + , 6_ respectively. We demon- 
strate numerically that for the case b + = 6_ and v + = u_ 
which, with appropriate cloice of the y — axis, repre- 
sents an odd function, 9 = 1/2, not 1/4 and we argue 
that 9 = 1/2 whenever there is no net drift, i.e. when 
t> + 6 + = V—b-. In these cases we argue that an effective 
diffusion in the x-direction, resulting from the stochas- 
tic movement between the alternating bands, dominates 
over the drift and naturally accounts for the exponent 
1/2. We obtain an analytical result in the limit 6_ — > 0, 
t>_ — > oo, with held fixed and equal to b+v+, and 

we verify that 9 = 1/2 in this case. 
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II. MODEL AND CALCULATION 



where 



The model we consider initially is defined by Eqs. (1) 
and (2), with the function f(y) given by 



«+ 



v+y 
-v-(-y) c 



y>o, 

y<0. 



(4) 



We work in the regime a + > a_ , where we expect a 
non-zero survival probability. The probability, Q(x,y) 
that the particle survives to infinite time satisfies the 
backward Fokkcr-Planck equation (BFPE) Dd yy Q + 
f(y)d x Q = 0. In the present context it is convenient 
to work instead with the 'killing probability', P = 1 — Q, 
which obviously satisfies the same equation, 



Dd yv P + f(y) d x P = Q, 



(5) 



but with different boundary conditions. The boundary 
conditions on P(x, y) are 



P(oo,j/) = 0, 
P(x,oo) = 0, 
P(0,y) = l, y<0, 



(6) 
(7) 
(8) 



the last of these following from the fact that, for an initial 
condition with x = and y < 0, the flow immediately 
takes the particle on to the absorbing boundary. 

Inserting the form (4) for f(y) in (5), the equation can 
be solved separately in each regime by separation of vari- 
ables. Imposing the boundary conditions (6) and (7) the 
result takes the form (setting D — 1 for convenience) 



P(x,y) 



,1/2/3-, 



J dka{k)K p+ {213+^k^y 1 

for y > 0, (9) 
V^J^dk {b(k)J _ (2/?_Vfc^(-y) 1 / 2 ' 3 -) 

+ c{k)J-p_ (2/3_ v ^T(-y) 1 / 2 ' 3 -)} e- kx , 

for y < 0, (10) 



where 



/3± = 1/(2 + a±), 



(11) 



J v {z) is a Bessel function and K v {z) a modified Bessel 
function. 

Eqs. (9) and (10) contain the three undetermined func- 
tions a(k), b(k) and c(k). The last two can be expressed 
in terms of a(k) by imposing suitable continuity condi- 
tions at y = 0. From the BFPE, Eq. (5), it follow that 
both P(x,y) and d y P(x,y) are continuous at y = 0. Im- 
posing these conditions gives the relations 



b{k) =K b k il3 +- l3 -^ 2 a{k), 
c(k)=K c k^--f 3 +V 2 a(k), 



(12) 
(13) 



K h = 



r(i + /3_) 



13+ 



2sin(/3 +7 r) r(l + /?+) (/J.^r)' 



2 r (/3+)r(i 



/?-) 



p+ ■ 



(14) 



(15) 



Finally, a{k) is determined by the boundary condition 
(8). This gives 



1 = 



y/=y J" dk [b(k)J p _ (2/3_V^(-J/) 1/2 ^) 

+ c{k)j_ (j _ (2/?_v^(-y) 1/2/3 -)} , V < 0- (16) 



In principle, Eqs. (12), (13) and (16) completely deter- 
mine the functions a(k), b(k) and c(k). In practice, how- 
ever, we have been unable to extract these functions ex- 
plicitly. We can, however, determine the small-/c forms. 
These are sufficient to demonstrate our main claims. 

To determine the small-A: behavior of a(k) we consider 
Eq. (16) in the limit y — > — oo, since the integral is dom- 
inated by small values of k in that limit. We make the 
ansatz 



a(k) ~ Ak a , k -> . 



(17) 



Then the small-A; forms of b(k) and c(k) are fixed by 
equations (12) and (13). Since /?+ < /L, 6(fc) > c(k) 
for fe — > 0, so the first term, involving b(k), in Eq. (16) 
dominates for large negative y, and the term involving 
c(k) is negligible. In this limit, therefore, Eq. (16) reads 
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AK b lim 

y— >-cx 



/ dkk a+ ^-f 3 -^ 2 

Jo 

xj f3 _ (2/?_V^(-y) 1/2 ^)- (is) 

Evaluating the integral fixes the values of the exponent 
a and the amplitude A in Eq. (17): 



A = 



(/?- 



(19) 
(20) 



Inserting the small-A; results for a(A), b(k) and c(fc) in 
Eqs. (9) and (10) determines P(x,y) in the entire regime 
where x or y (or both) are large. Especially simple is the 
limit of large x at y = 0, where one obtains 



P(x,0) 



r(/3+)r(/3_ - /3+) 03- v^) 



/3- 



2^ 6 r(/?_) (/? +v ^+ 



(21) 



for x — > oo. This result demonstrates that there is a non- 
zero infinite-time survival probability provided /3_ > /3+ , 
i.e. a+ > a_. The spatial decay exponent, /3_ — 
should have (we feel) a simple physical derivation, but so 
far it has eluded us. 
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III. PERIODIC FLOW FIELDS 

In the remaining part of this paper we consider the 
case where f(y) is a periodic function of y 7 with period 
2b. In particular, we consider the case 



v, 2nb < y < (2n + 1)6, 
—v, (2n - 1)6 < y < 2nb. 



(22) 



for all integers n. This function f(y) is odd, but we shall 
show that the survival probability decays asymptotically 
ast- 1 /?. 

We first present, in Figure 1, the result of a numeri- 
cal simulation for different values of the width, 6, of the 
alternating bands. In each case the particle is initially 
located at x = 0, on the interface between two bands. 
To avoid immediate absorption, the first step is taken 
into the band with positive velocity. For early times, 
t <C b 2 /D, the particle explores the two bands either side 
of the initial position and the decay exponent 1/4, ap- 
propriate to infinitely wide bands is obtained. At time 
t ~ b 2 /D, however, where the particle starts to explore 
many bands, there is a clear crossover from the initial 
i -1 / 4 decay to an asymptotic i -1 / 2 decay. 



-10 




10 
ln(t) 



15 



FIG. 1. Double logarithmic plot of the survival probabil- 
ity Q against time t for a range of band widths b. The data 
show, after a short-time transient, an initial £ _1//4 decay cross- 
ing over to a t' 1 ^ 2 decay at later times. The dashed lines with 
slopes —1/4 and —1/2 are guides to the eye. 

To understand these results, we first recall the use of 
the Sparre Anderson theorem [6] to demonstrate that 
9 = 1/4 for odd functions f(y). First suppose that 
f(y) > for y > and f(y) < for y < 0. We focus 
on crossings of the x-axis and regard the steps between 
crossings as the elementary steps of a one-dimensional 
random walk on the x-axis. The steps are clearly (i) in- 
dependent, and (ii) drawn from a symmetric distribution. 
The Sparre Anderson theorem states that the probability 



that the process has not crossed the absorbing boundary 
at x = after AT steps decays as Qn ~ A -1 / 2 for large 
A. Since the number of crossings in time t scales as 
A - i 1 / 2 , this implies Q(t) ~ i" 1 / 4 . 

Let us examine this argument more carefully. It's va- 
lidity clearly requires that a crossing of the absorbing 
boundary at x = can occur at most once between two 
consecutive crossings of the x-axis. This will always be 
the case if f(y) takes only one sign for y > and the 
other sign for y < 0. In the periodic case, however, f(y) 
takes both signs for y > and both signs for y < 0. The 
process can therefore cross the absorbing boundary, from 
x > to x < 0, and cross back again without crossing 
the x-axis. Such absorption process will be missed in a 
description which focuses only on crossings of the x-axis. 
Such a description will typically underestimate the value 
of 9 (it gives a lower bound on 9). 

A better way of looking at this problem is as follows. 
The probability Q(x,y,t) that the particle survives until 
time t, given that it started at (x, y), is periodic with pe- 
riod 26. In particular, the partial derivative d y Q vanishes 
by symmetry at the center of any band. It therefore suf- 
fices to solve the problem in the region —6/2 < y < 6/2, 
x > 0, which defines a semi-infinite strip, with boundary 
conditions Q(0,y,t) = 0, d y Q(x,±b/2,t) = 0. The flow 
velocity is v in the upper part of the strip (y > 0) and —v 
in the lower part. This strip problem has been discussed 
by Redner and Krapivsky [2], who argue as follows. The 
typical time between crossings of the center line, y = 0, 
is r ~ b 2 /D. The typical distance travelled, in the x- 
direction, in this time is I ~ vb 2 /D. The two cases x 3> I 
and x -C I, where x is the initial displacement of the 
particle, have to be analysed separately. The following 
discussion is based on reference [2] . 



A. The regime x 2> I 

The crossings of the x-axis define a symmetric one- 
dimensional random walk with step length ~ I and time 
step ~ r. The effective diffusion constant for motion in 
the x-direction is, therefore, Dm ~ 1 2 /t ~ v 2 b 2 /D. The 
survival probability for this random walk is given by the 
well-known result [4] 



Q(x,y,t) 




(23) 



Note that the y-dependence is not determined by this ar- 
gument. It is intuitively clear, however, that the result 
is insensitive to y at large x. 

Comparing this result with Figure 1, we note that the 
amplitude of the i -1 / 2 behaviour in the late-time regime 
increases with increasing 6 in Figure 1, whereas Eq. (23) 
predicts a decreasing amplitude with increasing 6. The 
data, however, were obtained with the initial x equal to 
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zero (with x = v after one time step). To make mean- 
ingful comparisons with the data, therefore, we need to 
analyse the opposite regime x <C I. 

B. The regime x -C I 

For x <C I, one can argue as follows. At short 
times, t -C t, the behaviour will be the same as for 
strips of infinite width, Q(t) ~ (ie/i^) 1 / 4 [2]. At later 
times, this result will be modified by a function of t/r, 
Q(t) = (x/vtf^fiDt/b 2 ). For t > r = fo 2 /L», there 
should be a crossover to the i -1 / 2 decay calculated above, 
since if the particle survives for time r, it will typically 
reach a distance x ~ vt = I. The scaling function f(z) 
therefore behaves as z -1 / 2 for large z, giving 

/xb 2 \ 1/4 1 

Q( *' M) ~Ud) ^' *» 52 /^- (24) 

The two results Eqs. (23) and (24) match, as required, 
at x ~ I = vb 2 /D, where both reduce to Q ~ b/y/Dt. 
The scaling form Q(t) ~ t' 1 ^ f(Dt/b 2 ) is tested in Fig- 
ure 2, where t x I^Q is plotted against t/b 2 on a double 
logarithmic plot. After the initial transients, there is a 
good collapse of the data onto a scaling curve, in which 
an initial slope of zero gives way, at t ~ b 2 , to a final slope 
of —1/4, corresponding to the regimes Q(t) ~ t- 1 / 4 and 
Q(t) ~ i -1 / 2 respectively. 

o | . 1 . 1 . 1 . 1 




-4 1 1 1 ' 1 1 1 ' 1 

-10 -5 5 10 

ln(t/b 2 ) 

FIG. 2. Double logarithmic plot of the survival probability 
Q against time t plotted in the scaling form £ 1//4 Q against 
t/b 2 . The dashed line with slopes —1/4 is a guide to the eye. 

We can interpret the asymptotic result Q ~ £ -1 / 2 as 
follows. If we coarse grain in the y-direction, the average 
deterministic drift is zero, and the nominally subdomi- 
nant diffusive motion in the x-direction [2], with diffu- 
sion constant Dm , plays an important role. The effective 



coarse-grained dynamics corresponds to anisotropic two- 
dimensional diffusion. If the absorbing boundary is the 
whole y-axis, diffusion in the y-direction is irrelevant (ex- 
cept for its role in inducing diffusion in the x-direction) 
and the standard one-dimensional result for the survival 
probability, Eq. (23), is obtained. One can now discuss, 
however, other types of absorbing boundary. First we 
consider absorption on the half-line x = 0, y < 0. Numer- 
ical simulations, presented in Figure 3, suggest 9=1/4 
for this case, while for absorption on the two half lines 
x = 0, y < and y = 0, x < 0, which define an absorbing 
wedge of angle 7r/2, the data suggest (for x > initially) 
61=1/3. 

These seemingly mysterious exponents are readily un- 
derstood in the context of the underlying anisotropic two- 
dimensional diffusion process. For isotropic diffusion in- 
side a wedge of opening angle <f>, with absorbing bound- 
aries on the edges of the wedge, the survival probability 
is known to decay as a power law, Q(t) ~ t~ e , with 
8 = Ti/2(j) [7,4]. The absorbing boundary on the half-line 
x = 0, y < corresponds to opening angle <f> = 27r, giv- 
ing 6 = 1/4, while the absorbing boundary on the two 
half- lines x — 0, y < and y = 0, x < corresponds 
to (j> = 37r/2, giving 9 = 1/3. Although the effective dif- 
fusion process is here anisotropic, it can be transformed 
into an isotropic process by rescaling the x or y variable. 
Since this rescaling does not change the geometry of the 
two absorbing boundaries considered, the results 9=1/4 
and 1/3 respectively are not changed. 



o 




_10 i . 1 . 1 . 1 . 1 

5 10 15 20 

met) 

FIG. 3. Double logarithmic plot of the survival probability 
Q against time t for the band model, for different types of ab- 
sorbing boundary. The dashed lines with slope —1/4, —1/3 
and —1/2 are guides to the eye. In order of increasing steep- 
ness, the corresponding data refer to absorbing boundaries 
given by (i) the negative j/-axis, (ii) the negative j/-axis and 
the negative a;-axis and (iii) the whole y-axis. Least squares 
fits to the data give slopes -0.254, -0.333 and -0.501 respec- 
tively, consistent with the predicted values. 
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It is straightforward to apply the same general reason- 
ing to the band model where the positive- and negative- 
flow bands are inequivalent. Consider a model where 
bands of width b + and f(y) = v + alternate with bands 
of width 6_ and f(y) = — u_. For the special cases where 
b + v + = 6_u_ there will be no macroscopic flow, and 
a coarse-grained description will lead to an anisotropic 
diffusion model as before, with the same values of 9, 
e.g. 9 = 1/2 if the whole y-axis is absorbing (indeed, 
9=1/2 will hold for any periodic f(y) with zero mean). 
If b + v + — 6_v_ > 0, the coarse-grained model has a 
mean flow away from the absorbing boundary, and the 
infinite-time survival probability is non-zero, while for 
b + v + — 6_u_ < the mean flow is towards the boundary 
and the survival probability decays exponentially with 
time. In the following subsection we provide an exact 
solution for the case b+v + = 6_v_, in the limit 6_ — > 0, 
V- — > oo, with u = 6_-u_ held fixed. 



C. A solvable model 

The model we study is defined by a periodic function 
f(y), with period b, given by 



f(y) 



>,6(y), -b/2<y<b/2, 



(25) 



with f(y + b) = f(y) for all y. The marginal case of inter- 
est, corresponding to no net flow, is given by u — bv. It 
is convenient, in the first instance, to compute the killing 
probability P(x,y,t). The periodicity implies that we 
need only consider the regime —6/2 < y < 6/2, with the 
boundary conditions d y P(x,b/2,t) = d y P(x, — 6/2, t) = 
0. The function P obeys the BFPE 



t P = Dd vv P +[v- u5(y)]d x P. 



(26) 



Taking the Laplace transform with respect to time, and 
exploiting the initial condition P(x,y, 0) = 1 for all 
x > 0, gives 



Dd yv P +[v- u5{y)]d x P -sP = 0, 



(27) 



where P(x,y,s) = J °° dt exp(—st)P(x,y,t). The gen- 
eral solution satisfying the given boundary conditions at 
y = ±6/2 is (setting D = 1 for convenience) 



P{x,y,s) 



f 

Jo 



dk a(k, s)e 
x cosh 



-kx 



Vkv + s i - 



(28) 



for i/^0. Imposing the discontinuity in d y P at y = 
implied by Eq. (27), d y P(x, 0+, s) - d y P(x,0-,s) = 
ud x P(x, 0, s) for all x > and all s, leads to the con- 
dition 



tanh7 = 



ub / 7 



7« \ b 2 



where 



7 = -Vkv + s. 



(29) 



(30) 



Eq. (29) shows that only a single value of k is possible 
for each s. Specialising to the marginal case u = bv, this 
equation simplifies to 



tanh7 



_ 6^s 

' 7 4 7 



(31) 



The large-time behavior is governed by the small-s so- 
lution of Eq. (31): 7 ~ (36 2 s/4) 1 / 4 . Inserting this in to 
Eq. (30) gives, for small s, k ~ 2^/3s/bv and \Jkv + s ~ 
{isf^^/tjb. Putting these results into Eq. (28) gives 

2v/3i \ 



P(x, y, s) = a(s) exp 



bv 



x cosh 



(3s) 



1/4 



6/2 



(32) 



for s — > 0. 

The amplitude a(s) is fixed by the boundary condi- 
tion, P(0,0,i) = 1 for all t, which gives P(0,0,s) = l/.s 
and a(s) = {s cosh[(3s) 1 / 4 y / 6/2]} _1 for small s. Since 
Q(x,y,t) = 1 - P(x,y,t), we have Q(x,y,s) = 1/s - 
P(x, y, s). The small-s result for P determines the large- 
t form of Q(x, y, t) as 



Q{x,y,t) = 



2V3 



bVnt 



\y\) 



(33) 



for \y\ < 6/2. This result confirms, for this special case, 
the result Q(t) <~ t^ 1 / 2 obtained above using general 
heuristic arguments for the case where there is no net 
flow. Note that Eq. (33) has the same form as Eq. (23), 
and that this form holds for all values of x/l, where 
/ = vb 2 / D, i.e. there is no analogue, in this model, of 
Eq. (24) for x » I. This is because in the present model 
there is no regime in which the systems resembles a two- 
band model. 



IV. DISCUSSION AND SUMMARY 

In this paper we have obtained further results on the 
first-passage properties of a particle diffusing in the y- 
direction and subject to a deterministic drift velocity, 
f(y), in the x-direction, with an absorbing boundary at 
x = 0. 

In the first part of the paper we considered the case 
where f(y) oc y a + for y > and f(y) cx — for 
y < 0. We confirmed the conjecture proposed in [1] that 
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the larger power is 'strongly dominating' (in contrast to 
the 'weakly dominating' behaviour obtained when the 
exponents a + and a_ are equal but the amplitudes are 
different [1]). We showed explicitly that for a + > a_ 
the survival probability approaches a non-zero limiting 
value, and we obtained an expression for this value. 

In the second part of the paper we considered the case 
where f(y) has a banded periodic form, with alternating 
bands of equal and opposite drift velocity. We presented 
numerical data, supported by a heuristic argument, to 
show that the asymptotic persistence exponent is 6 = 1/2 
in this case, and obtained a generalisation for different 
forms of absorbing boundary. We argued that the re- 
sult generalises to any model where the coarse-grained 
flow field vanishes, and exemplified this result through a 
soluble model. 



PCs work was supported by EPSRC. 



[1] A. J. Bray and P. Gonos, J. Phys. A 37, L361 (2004). 
[2] S. Redner and P. L. Krapivsky, J. Stat. Phys. 82, 999 
(1996). 

[3] T. W. Burkhardt, J. Phys. A 33, L429 (2000). 

[4] S. Redner, A Guide to First-Passage Processes (CUP, 

Cambridge, 2001). 
[5] S. N. Majumdar, Curr. Sci. India 77, 370 (1999). 
[6] E. Sparre Andersen, Math. Scand. 1, 263 (1953); ibid. 2, 

195 (1953). 

[7] M. E. Fisher and M. P. Gelfand, J. Stat. Phys. 53, 175 
(1988). 



G 



